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ABSTRACT 

We present a method for deriving a smoothed estimate of the peculiar velocity field 
of a set of galaxies with measured circular velocities r] = logAu and apparent magnitudes 
m. The method is based on minimizing the scatter of a linear inverse TuUy-Fisher relation 
r] = ri[M) where the absolute magnitude of each galaxy is inferred from its redshift z^ 
corrected by a peculiar velocity field, M oc m — blog[z — u). We describe the radial 
peculiar velocity field u[z) in terms of a set of orthogonal functions which can be derived 
from any convenient basis set; as an example we take them to be linear combinations of 
low order spherical harmonic and spherical Bessel functions. The model parameters are 
then found by maximizing the likelihood function for measuring a set of observed rj. The 
predicted peculiar velocities are free of Malmquist bias in the absence of multi-streaming, 
provided no selection criteria are imposed on the measurement of circular velocities. This 
procedure can be considered as a generalized smoothing algorithm of the peculiar velocity 
field, and is particularly useful for comparison to the smoothed gravity field derived from 
full-sky galaxy redshift catalogs such as the IRAS surveys. We demonstrate the technique 
using a catalog of "galaxies" derived from an N-body simulation. Increasing the resolution 
of the velocity smoothing beyond a certain level degrades the correlation of fitted velocities 
against the velocities calculated from linear theory methods, which have finite resolution, 
but the slope of the scatter diagram, and therefore the derived density parameter, remains 
fixed. 



Subject headings: cosmology — dark matter — galaxies: clustering — galaxies: formation 
— gravity — large-scale structure of the universe 



1. INTRODUCTION 



The large scale motions of galaxies in the universe are an interesting probe of the 
nature of large scale structure in the Universe. Indeed, if galaxies move under the influence 
of the gravity of the underlying matter fleld, then their peculiar velocities, the deviations 
from smooth Hubble expansion, can be used to address fundamental questions such as 
the amplitude of the clustering in the dark matter distribution, the mean density of the 
Universe, and the statistical nature of the primordial fluctuations. A thorough review of 
this active fleld is given by Dekel (1994). 

All of the methods of mapping deviations from Hubble flow depend on an observable 
parameter that is distance dependent, e.g. magnitude, and another which is distance 
independent, e.g. line width or rotational velocity. The remarkable TuUy-Fisher (TF) 
correlation between rotational velocity and absolute magnitude for disklike galaxies is the 
most commonly used distance indicator, with a distance error per galaxy of order 16%. 
Alternative extragalactic distance estimators have been recently reviewed by Jacoby et al. 
(1992). 

In a widely cited recent publication, Dekel et al. (1993) compared the density fleld of 
IRAS selected galaxies with the divergence of the measured peculiar velocity fleld derived 
from the POTENT algorithm (Bertschinger et al. 1990; Dekel, Bertschinger, and Faber 
1990). They have shown that the flows are directed toward the peaks of the galaxy density, 
as expected in the linear theory of gravitational instability (Peebles 1980); quantitatively, 
they derive a high estimate of the density parameter, j3 = Q^'^ /bj sa 1.3 it 0.3, where bj is 
the linear bias in the IRAS galaxy distribution relative to the mass distribution. 

The Dekel et al. (1993) result depended on a forward flt (i.e. M = M[ri))^ of the TuUy 
Fisher relation for spiral galaxies, and the equivalent D^ — a relationship for the elliptical 
galaxies, and is subject to homogeneous as well as inhomogeneous Malmquist bias (see 
e.g. Lynden-Bell et al. 1988, Gould 1993, Willick 1994). These effects must be carefully 
calibrated, and are quite sensitive to the estimated uncertainty in the TF scatter. The 
POTENT method requires extensive smoothing, as well as the gridding of the peculiar 
velocity data, which can introduce a serious "sampling gradient bias" . 

Such a complex analysis must be tested by a number of independent means, using the 
best available data. Willick et al. (1994b) are using the newly compiled Mark3 dataset 
(Willick et al. 1994a) of sa 2800 spiral galaxies with high quality TuUy-Fisher parameters 
(i.e. apparent magnitudes and rotational velocities), to compare the peculiar velocity pre- 
dicted from the distribution of IRAS galaxies to that observed for individual objects, using 
a maximum likelihood method that takes into consideration the effects of an inhomoge- 
neous mass distribution and a variety of sampling selection effects (Willick 1994). Nusser 
and Davis (1994) used an early version of this Mark3 dataset processed through the PO- 
TENT algorithm to compare the inferred dipole of the radial peculiar velocity fleld to that 
expected from the 1.2 Jy IRAS redshift survey data (Fisher et al. 1994a). They derived a 
best estimate fi sa 0.6 but they did not attempt to estimate the uncertainty because they 
had no control over the POTENT fleld. 



Direct estimation of the peculiar velocity of individual galaxies based on the forward 
TuUy-Fisher relationship can be used to generate a peculiar velocity field, with nice visual 
images, but is difficult to correct for the myriad selection biases. Less problematic in terms 
of Malmquist corrections is to use the inverse TuUy Fisher relationship (Schechter 1980; 
Aaronson et al. 1982; Lynden-Bell 1991), (i.e. rj = rji^M)), but this relationship is difficult 
to visualize, and is best done by means of a parameterized model. Aaronson et al. (1982) 
fitted a spherically symmetric Virgocentric flow model to their data using an inverse TF 
relationship and Lahav (1991) fltted a pure bulk flow model. Roth (1993) and Willick et 
al. (1994b) flt a series of IRAS predicted flow models (Strauss and Davis 1988; Yahil 1988) 
to their data where fi used to construct the IRAS flow models is the only free parameter. 
Although this method can derive a 'best' estimate of /3, it does not immediately tell us 
if the model is a good flt to the data. This procedure also suffers if the redshift-distance 
relation is multivalued, as frequently occurs toward clusters of galaxies. 

In this paper we shall demonstrate a middle ground, in which, without binning the 
data, we effectively generate a smoothed version of the peculiar velocity fleld by means 
of an inverse flt to the TF relationship. A smoothed flt to the observed data has the 
advantage that it fllters out small scale peculiar velocities that are beyond the resolution 
of the IRAS gravity maps, and which are not expected to obey linear theory in any event. 
The result is that, without ever binning the peculiar velocity data, we can generate maps 
of the peculiar velocity fleld, which we demonstrate are an unbiased estimate of the true 
fleld. 

In section 2, we describe the method, its biases, and the construction of the orthogonal 
functions. In section 3, we construct a Monte-Carlo catalog derived from an N-body 
simulation to test the method, and demonstrate that the predicted velocity fleld is unbiased 
for different choices of resolution. A discussion of the utility of the technique is given in 
section 4. We shall apply this method to the Mark3 data in a separate publication. 

2. THE METHOD 

2.1 Maximum Likelihood fit for r] 

The goal is to derive a spatially smoothed estimate of the peculiar velocities of a 
sample of galaxies having measured circular velocity parameters, rji^ apparent magnitudes, 
nii^ and redshifts, Zi. We assume that the circular velocity parameter, ry, of a galaxy is 
related to its absolute magnitude, M, by means of a linear inverse TuUy-Fisher (ITF) 
relation, i.e., 

r] = sM + r]o+R, (1) 

where rjo and s are constants and i? is a Gaussian random component with mean and 
rms value cr^ independent of r] and M. If u is the radial peculiar velocity of a galaxy then 
the apparent and absolute magnitudes are related by, 

M = m- 5log{z -u) - 15. (2) 

where the peculiar velocity and the redshift have units of km/s. The relations (1) and (2) 
suggest the following procedure for estimating the peculiar velocities: 



i) Express the peculiar velocities in terms of a parameterized model. 

ii) Use the peculiar velocity model to infer absolute magnitudes, Mi (ttz^, 2;^, model), 
of galaxies from their apparent magnitudes and redshifts given the relation (2). 

iii) Assume Gaussian deviations from the ITF relation so that the probability of 
measuring the set of observed values rji^ given (1) and the model predicted absolute 
magnitudes, M^, given from ( ii), is: 



Pr(r7|M) oc exp 



rio + sMj - r]A 



2 



U 2..^ 



where Ng is the number of galaxies in the sample. 

( iv) Find the model best fit parameters by maximizing the probability function (3). 

It should be noted that a peculiar velocity model in which the number of fit parameters 
is significantly less than the number of galaxies in the sample is fully consistent with the 
gravitational instability process, where large scale, coherent flows are expected to have 
substantial amplitude. 

The procedure for flnding the best flt parameters for the model is equivalent to mini- 
mizing the x^ function, 

(r/o + sM^ - r]A 



"i 



where cr^ is constant. We now proceed to construct a form for the flt model. Using (2) we 
write, 

M, = Mo. + P. (5) 



where 



Moi = m — 5log[zi) — 15 and Pi = —Slog (1 ^ ) . (6) 



Since the absolute magnitudes depend linearly on the quantities P^, we choose to express 
Pi in terms of linear combination of a set of independent functions. It is the parameters of 
this combination which we shall determine by minimizing the x^ . Thus we write for each 
galaxy. 



Jmax 

P^=Y1 «'^/ (7) 

J=0 



where jmax is the total number of the flt parameters a^ and the matrix F satisfles the 
following orthonormality condition. 



J^F^Ff=5i-\ 



i=l 



where 8k is the Kronecker deha and we set i^° = l/yNg. After substituting (5), (6) and 
(7) in (4) we find, 

Thanks to the orthogonahty condition (8) there is no couphng between the diff"erent a? 
modes. The TF parameters rjQ and s depend on the velocity model, therefore the x^ 
should be minimized subject to variations in s and r/o as well as in a? . At the minimum, 
the derivatives dx^ /da^ , dx^ /drjo and dx^ /ds vanish, yielding respectively, 

Y,{r]o + sMo^ - r],)Fi + sd.^ = 0, 

i 

J2irio + sMo^ -r]^) + sJ2 o' Ff = 0, 
i hi 



(10) 



hJ 



The first two of Eqs (10) are degenerate for j = 0. Then the determination of a° is 
somewhat arbitrary. The reason for this ambiguity is that (P ^ amounts to a radial 
peculiar velocity u oc z which is indistinguishable from the Hubble expansion. For a 
realistic database with finite surveyed volume about us, a non vanishing a° implies that 
the surveyed volume is a region of constant density contrast. This is highly improbable 
for most current observations which typically extend out to distances larger 8000 km/s. 
Therefore we set a° = 0. With this choice for the zeroth mode, the solution for a^ , rjo and 
s can be written in a closed form as, 

S-'=-^ J](3Mo. -r/,)^/, 

i 

^o=-^i:(3Mo.-,.), ^^^^ 

E,, Mo.ry, (l + iV, E, F^F^) - N, E. Mo.ry, 



Z^.Mo.Mo,(l + N,Z,K'H)-N,Z^M^,, 



The terms involving the matrix F are the correction in the estimate of the ITF parameter 
s arising from the peculiar velocities. The covariance matrix can be computed rigorously. 
However if the correction terms are small, we may neglect the coupling between a^ and s. 
There is no coupling between a^ and rjo because of the orthogonality of the modes and our 
choice of the zeroth mode. This implies that for Gaussian scatter of r] around the ITF, the 



errors are normal and given by, 

1 d\^ 




C^ai 



^m = I oTT-T ) = -^T^^ fl2) 



CTc 



It is not surprising that the errors in a^ scale like 1/s. Indeed, for s = the ITF is 
superfluous and cannot constrain the peculiar velocity fleld. Using the linear combination 
(7) we flnd that the error in the determination of the quantity Pi of any individual galaxy 

is, 

/. \ 1/2 

/ Jmax \ 

-'■. = ^\E(Fir\ (13) 



j=i 

The orthonormality condition (8) implies that ap^ oc ^/jmax/Ng. The corresponding error 
in the velocity in the limit u/ z <^ 1 is, 

(Ju,=0A6(jp,z, (14) 

In the limit jmax -^ Ng we recover the expected error, cr^; sa (jj^Zi/s. 

2.2 Biases in The Method 

The ITF method described above is intended to provide a generalized smoothed flt 
to the peculiar velocity data, and is not in itself motivated by any particular assumptions 
about the nature of the flow. We do not, for example, require the flow to be irrotational, 
as is required for POTENT. The flt makes no assumptions about the conversion of redshift 
space to real space, and simply seeks the best smoothed flt in redshift space. An important 
question to ask of any method is whether it suffers from any systematic biases, and if so, 
can they be readily calibrated and corrected. 

There are three effects which bias the estimation of peculiar velocities from the ITF. 
First, there might be some selection effects associated with measuring circular velocities. 
An example of this is the reduced signal to noise ratio of broad 21cm lines, relative to 
narrow lines of the same equivalent width. This is not considered to be a problem in 
practice, as galaxies with larger rotational width usually contain more hydrogen than 
those with slower rotation. Furthermore, many of the recent observations of line width are 
based on Ha rotation curves derived from long slit spectra, for which there is no signal to 
noise penalty for broad lines. 

A second, more important effect, is that small scale velocity dispersion introduces 
homogeneous as well as inhomogeneous Malmquist bias (e.g. Lynden-Bell 1991, Willick 

6 



1994). To see the effect of velocity dispersion, consider the following example. Suppose 
that we observe a set of [rji^mi) for Ng galaxies, all at the same redshift, z^ and direction 
in the sky. Assume that the parameters rjo and s of the ITF have already been determined. 
Let us for the sake of simplicity consider a case in which all the galaxies have a common 
radial peculiar velocity component, uncorrelated random components, 5ui^ with mean. 
We seek an estimate, Uest, for Utrue- By maximizing the likelihood for measuring the set 
of rji with respect to the estimate, Uest, we find, 

5sNg log (l - ^) = J](r7o + sMo^ - r,,) (15) 

i=l 

But Moi = Mi + 51og(l — Ui/z)^ where Ui = Utme + ^Ui are the peculiar velocities of the 
galaxies. Since ^j(?7o + -sMi — rji) = 0, we obtain, 

iV,log(l-^)=5]log(l-^) (16) 

i 

In the absence of small scale peculiar velocities, we have 5ui = and Ui = Utrue- Thus the 
last equation yields an unbiased Uest = Utrue- However, if 5ui ^ 0, then, by expanding the 
logarithmic terms to first order, we find. 



1 ^^ 



u 



est 



Utrue + ^2^Su, (17) 



'a 



i=l 



The sum of 5ui evaluated at a given redshift does not necessarily vanish on the average, 
and the resulting bias Aubias = ^ 5ui/Ng caused by the velocity dispersion can be written 
in the general form, 

J 5uPT{5u)PT{r = z — u)du 



Aubto 



J Pr {5u) Pr(r = z — u)d 



where we assume that the probability distribution of 5u is independent of position and that 
Pr(r), the probability that an object lies at a position r, is r'^n{r)^{r) where n{r) is the 
underlying number density and $(r) represents the selection probability that an object is 
actually identified in the sample. If 5u is normally distributed about Utrue with dispersion 
(T„, then, by expanding Pr(r = z — u) to first order in u about z^ Eq(18) yields, 

^ dlnl z'^n(z)^(z)] , , 

Aub^as = ol ^ . • 19 

zqXvlz 

Approximating the selection function in a finite region by a power law $ oc r~", we find, 

/ dlnn\ cr,^ 

Au6^a. = 2-a + -^^. 20 

\ din^; / z 

For typical values, Ou = 200 km/s, z = 2000 km/s, the bias is 40 km/s for a homogeneous 
distribution of objects in a volume limited catalog (n = const ^ a = 0). An inhomogeneous 
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spatial distribution of galaxies, however, would be expected to lead to larger bias in the 
components of the fitted flow. 

We point out that this same bias occurs with a larger effective value of cr^ in multi- 
valued zones, where objects with very different distances have similar redshifts, due to 
infall peculiar velocities toward a high density peak. Again, the bias is expected to depend 
inversely on the redshift. However the dependence on the density run and the selec- 
tion function is more diflicult to evaluate as the distribution of peculiar velocities in a 
multi-valued zone depends on the proflle of the density peak. In realistic applications, an 
assessment of this bias in the ITF estimated peculiar velocities may be made by inspecting 
the multi-valued zones seen in the prediction of velocity flelds from galaxy redshift surveys 
such as the 1.2Jy IRAS. 

In our methodology, the smoothed velocity fleld is deduced from the modeled function 
P. This is desirable since it avoids a possible bias in the calibration of the ITF parameters 
within magnitude limited surveys due to small scale velocity dispersion. Indeed, in mag- 
nitude limited surveys, objects at larger redshifts tend to have larger absolute luminosity. 
To second order in u/z, P sa 2.n{ul z -\- [u/zY /2). If the velocity fleld at redshift z is 
u, then averaging over a sample of galaxies with small scale dispersion cr^ at this position 
gives an average P of 

- ^ /" u u^ a^, 

P^2.17 - + — + -^ 

\z Zz^ Zz^ 

The flrst two terms represent the estimate of P in terms of the linear modeling of the 
peculiar velocity fleld, but the third term amounts to a bias that diminishes with distance, 
i.e. with absolute luminosity. Therefore the inferred slope of the ITF relation is expected 
to be biased by the presence of small scale velocity noise when expressing velocity fleld in 
a linear model. The bias would result in a flatter slope s if u is the true mean peculiar 
velocity. Since the model velocity fleld and ITF parameters are found self-consistently, 
this bias would affect the estimated peculiar velocities. However, by using a linear model 
for P rather than u, we avoid this form of calibration bias, but at the same time, a third 
bias of our estimate of the velocity is introduced. Our estimate of the peculiar velocity at 
a redshift z can, to second order in P, be written as, 

^^ = 0.46(P-0.23P^) . 

z 

But the true mean peculiar velocity is obtained by averaging the quantities Pi of all objects 
at the redshift z. Therefore the true mean peculiar velocity is, 

^^^^^ = 0.46 (P - 0.23P^ - 0.23a?>) , 

z 

where ap is the rms dispersion of P about P due to small scale velocity dispersion. And, 
for simplicity, we have ignored the bias in P caused by averaging at the redshift and not 
the true distance of objects. Hence our estimate of Uest is biased by the presence of the 
term with ap. To the relevant order, ap oc Uu/ z, hence the bias is oc cr'^/ z which is also the 



behavior of the bias (20). At large redshifts this bias is not important since there u oc P 
and u oc P. In our method the slope of the ITF as inferred from the solution (11) differs 
from the true one only because of the bias given by (20). 

Clearly, the biases associated with the small scale velocity dispersion discussed here 
can be avoided by restricting the application of the method to large enough redshifts. 

2.3 Construction of The Orthonornial Matrix 

To implement the method we need to construct the matrix F. In principle this can 
be done by performing a Gramm-Schmidt algorithm on any set of independent functions. 
In the case of a uniform distribution of matter, an orthogonal base of functions is formed 
from the radial spherical Bessel and the angular spherical harmonic functions (Regos & 
Szalay 1989, Fisher et al. 1994b). Here we shall use Bessel functions to describe the radial 
behavior of the angular multipoles obtained by the expansion in spherical harmonics. We 
write the quantity P in the form, 

f^max ''max ^^ — ^ 

P{z,e,4>)= Y.Y.Y. "-f^iJiiknz)-cn)Yi^ieA) (21) 

n=0 1=0 m= — l 

where jj is the derivative of the spherical Bessel of order /, Yim are the angular spherical 
harmonics, and anim is the amplitude of the mode with /, m and n. The factor cn is 
non-zero only for 1 = 1. The actual value of cn is determined by the reference frame 
relative to which the redshifts are measured and hence the peculiar velocities are estimated. 
The natural choice is to work in a frame in which the observer is at rest relative to the 
flow of matter near the origin (e.g. Nusser & Davis 1994). In the real universe this 
would be the Local Group frame. In such a frame, P vanishes near the origin, hence 
cii = linix^o j[{x). In the limit \u/z\ <^ 1, which is expected both near the origin and at 
the outer limit, we have P oc u/ z. Therefore the expansion (21) is well motivated because 
it corresponds to the correct asymptotic behavior at the origin of a velocity potential $, 
$ = Y.aimnji{kn)Yim{6,(l)) (Jacksou 1962). 

Therefore we set the basis set of functions to be orthonormalized on the discrete 
non-uniform distribution of the data points as, 

i^/ = - [j'l {knz^) - cn) Yi^ (^„ </>,) , (22) 

z 

where j = l...jmax is deflned by each triplet n,/, and m. It is those functions which 
we orthonormalize using the Gramm-Schmidt algorithm. The orthonormal functions are 
denoted by F. 

The expansion of an arbitrary radial function into a set of Bessel functions requires the 
speciflcation of additional boundary conditions appropriate for the problem in question. 
These boundary conditions determine the wavenumbers k^ in the expansion (21). In any 
realistic application, it is necessary to set a maximum radius, z^ax^ within which the data 
is sufliciently reliable to constrain the peculiar velocity fleld. It is at z^ax where the outer 



boundary conditions are specified. We are guided by the asymptotic behavior of a pecuhar 
velocity field generated by fluctuations in the mass-density distribution. If we assume that 
the distribution of matter beyond z^ax is uniform, then potential theory dictates that 
u oc 2;~('+2) Sit z > Zjnax- In addition to a decaying mode z~'^ , the dipole harmonic, / = 1, 
attains a constant value which is the reflection of the observer's motion relative to the bulk 
motion of matter at z > z^ax- Thus we require 

dhl[jl[knZmax) _ // , r,\ (^n\ 

dln(._) - ^^ + ^^' ^^^^ 

which, using the recurrence relations of spherical Bessel functions, leads to ji-i{knZmax) = 
(Fisher et al. 1994b). We solve for the lowest values of kn which satisfy this relation for 
each /. 

Once the zeros are flxed we introduce additional quadrupole modes which take into 
account the velocity components generated at z < z^ax by mass fluctuations at z > z^ax- 
An external quadrupole fleld generates a peculiar velocity which scales as z; External fluc- 
tuations do not generate dipole and monopole velocity terms at z < z^ax^ and we neglect 
higher order modes. This external quadrupole mode is not orthogonal to the oscillatory 
Bessel function modes, but this does not prevent including it within a complete basis set 
of functions. The external source terms are expected to make a negligible contribution to 
the predicted velocity fleld in the internal region unless the external density fluctuations 
have a gross quadrupole term. 

The total number of the model flt parameters is flxed by the highest order of the 
spherical harmonic, Imax and the number of radial modes in the Bessel function expansion, 
Umax- Including the external modes for / = 2, we have 

^modes — ^rnax\''max ~r '^'"max ~r J-J ~r V^^J 

valid for Imax ^ 2. Thus for a typical case, rimax = 8 and Imax = 3 we have 133 degrees of 
freedom in the flt. 

3. TESTING WITH AN N-BODY SIMULATION 

We test the reconstruction of peculiar velocities from the ITF using a standard CDM 
N-body simulation normalized to have rms mass fluctuations — = 0.7 in a sphere of 
8h~"^Mpc. Particle velocities in the simulation were smoothed with a Gaussian window 
of width lh~"^Mpc so that the pair velocity dispersion on a scale of lh~"^Mpc is close 
to the observed value, 317km/s (Fisher et al. 1994a). A particle moving at 600 km/s, 
with small local shear, and not in a dense region, was chosen as the central observer 
(Fisher et al. 1994a). The coordinates were rotated so that this observer is moving in the 
same direction of the sky as our Local Group of Galaxies, according to the CMBR dipole 
anisotropy. A mock full sky 1.2Jy IRAS catalog was generated from the simulation by 
assigning luminosities to particles according to the luminosity function of the 1.2Jy IRAS 
survey. This procedure resulted in 4881 "galaxies" out to a distance of ISOh""^ from the 
"observer" . 
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To further improve the reahsm of the simulation, we scanned the galaxy positions of 
the Marks catalog (Willick et al. 1994a), which is quite anisotropically selected across the 
sky, and selected the nearest simulation "galaxy" to each MarkS real galaxy, resulting in 
a final catalog of 2871 "galaxies". An inverse TF relation with r/o = —2 and s = —0.1 was 
used to generate three mock MarkS catalogs with cr^ = 0.01, 0.03 and 0.07. 

In Figure 1 we show the reconstruction of the peculiar velocities from the ITF relation 
using the fake circular velocities and apparent magnitudes of the 4881 "galaxies" in the 
mock IRAS catalog with cr^ = 0.01. The reconstruction method was applied once with 
{iT-max} ^max) = (8, 3) and ouce with {rimax} ^max) = (12, 4) radial and angular modes respec- 
tively. The maximum redshift at which the external boundary conditions used to determine 
the zeros of the Bessel functions was z^ax = 10000 km/s. We see that the method is suc- 
cessful at recovering the true velocities of "galaxies" in the simulation. The reconstruction 
with (nmaxT^max) = (12,4) shows more scatter than the one with (nmaxT^max) = (8,3). 
The reason is that the errors in the estimated velocities increase with the total number of 
fit parameters. In other words, though the x^ improves as the number of fit parameters is 
increased, it attains a wider minimum. The difference between the estimated and the true 
peculiar velocities as a function of redshift is shown in Figure 2. As expected the scatter 
increases with redshift. Yet, no biases are detected in the estimated versus true velocities 
at any redshift. The predicted slope of the ITF relation as given from the solution (11) is 
s = —0.098 while the true value is s = —0.1, demonstrating that this bias is negligible. 

Most of the scatter seen in Figures 1 and 2 is a result of the small scale velocity 
dispersion which is not fitted by the smoothed field, since the scatter in the ITF, cr^ = 0.01, 
translates to a magnitude error of (Tm = cr^/s = 0.1, or a peculiar velocity error of 0.05^;, 
considerably smaller than is seen in Figure 2. 

Figure 3 shows the reconstruction of the true peculiar velocities of the mock Mark3 
galaxies for the lower resolution smoothing. Only velocities of "galaxies" within a redshift 
of 5000 km/s are plotted. The two upper and the lower right panels of Figure 3 show 
the reconstruction of the true peculiar velocities of the mock Mark3 data for the three 
values of cr^. Most encouraging is that no systematic bias in the estimate of the peculiar 
velocities is evident. The scatter in the plots for cr^ = 0.03 and 0.07 is attributed to the 
scatter in the ITF relation and the small scale intrinsic velocity dispersion. The predicted 
slopes of the ITF are s = -0.0978, -0.0974 and -0.0964 for a^ = 0.01, 0.03 and 0.07. To 
measure how well our smoothed estimate of the peculiar velocities can match the velocity 
field as predicted from the redshift space distribution of galaxies, we show in the lower left 
panel the smoothed peculiar velocities with cr^ = 0.03 against the linear theory velocities 
predicted from the density field in the mock IRAS catalog. The velocities from the density 
field were obtained using the linear theory method of Nusser & Davis (1994). It seems 
that the reconstruction of the velocities from the simulated ITF and the density field agree 
very well, with no apparent bias in the slope of the correlation. This result motivates the 
application of this smoothing method to real data so as to better constrain the value of 
f3 = 0°'^/6. The difference between the Y and X axes of each panel in Figure 3 is shown 
versus redshift in Figure 4. Again we note that even in sparse data like the Mark3 catalog, 
the method is capable of recovering the peculiar velocities without any noticeable biases. 
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Figure 5 shows the ratio of the velocity differences from Figure 4 to the error in the 
estimated velocities versus redshift. The errors were approximated using Eq(13). Since 
(Tu oc z^ the scatter decreases with redshift. In the absence of any small scale velocity 
dispersion, the rms value of the scatter in this plot would be unity. The large scatter 
for cr^ = 0.01 indicates that small scale velocities amounts to most of the random error 
in the estimate of peculiar velocities of individual "galaxies" . But the scatter caused by 
small scale velocity dispersion is not negligible even for cr^ as large as 0.07. No systematic 
bias of the estimated relative to the true peculiar velocities is evident. However, the ratio 
corresponding to the linear theory prediction plotted in the lower left panel shows some 
positive offset. 

In Figure 6a and 6b we examine in more detail the nature of the agreement between the 
smoothed and the true peculiar velocities of the mock MarkS "galaxies", with cr^ = 0.03. 
We show a full sky plot of the positions of all galaxies with 3500 < z < 5000 km/s. 
The differences between the estimated and true velocities are in very good agreement 
considering that the estimated velocities are smoothed on a scale ~ 1000 km/s. 

In order to demonstrate the possible biases which might be present in estimated 
peculiar velocities, we have boosted the small scale velocity dispersion by assigning to 
"galaxies" in the mock IRAS catalog additional radial peculiar velocity components drawn 
from a normal distribution with dispersion cr^ = 400km/s so that the total small scale 
velocity dispersion is ~ 500km/s. The new peculiar velocities were then used to compute 
the redshifts of "galaxies" . In Figure 7 we plot the difference between the estimated and 
true peculiar velocities versus redshift. We used cr^ = 0.01 in this reconstruction. It is 
clear the the estimated velocities are biased. But even for such a large value of the velocity 
dispersion the bias is not severe. The sense of the bias is that at large and low redshifts 
(z > 4000km/s, z < 2000km/s) the method underestimates the true peculiar velocities. 
The bias at large redshifts is caused by the lower slope predicted by the method. The 
predicted slope using the solution (11) is s = —0.092 while the true value is s = —0.1. 

As an indication of the coherence of the residuals of the fit Au^, we display in Figure 8 
the correlation function of the velocity residuals, computed in the same fashion as Gorski 
et al. (1989). We plot the quantity 

y;,- i W^W./^u, Au,cos(6'„) 

where the sum is over all pairs in the mock Mark3 catalog (cr^ = 0.03) separated by 
redshift separation As, cos[6ij) is the angle on the sky between the galaxies i and j, and 
the weight function is arbitrary but is here chosen to be unity, Wi = 1 for z < 6000 km/s 
and Wi = for z > 6000 km/s. The short dashed curve in Figure 8 shows |\I'i,p^ for 
the true peculiar velocity of the mock Mark3 data, not the residuals; in the local group 
frame shown here, the field is almost completely coherent. The solid line shows |\I'i,p^ 
(with sign preserved) of the fitted minus true peculiar velocity for the Mark3 mock sample, 
while the dashed curve shows |\I'i,p^ for the fitted minus linear theory peculiar velocity. 
Note in the latter case that the coherence of the residuals is considerably larger than for 
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the fit-true velocities. Tliis colierence is consistent with the small offset of residuals from 
zero seen in the lower left panel of figure 5 (compare to upper right panel of figure 5) and 
is due to the nonlocal nature of the velocity field inferred from linear theory. Errors in 
this field due to shot noise will result in coherent velocity residuals. The extra coherence 
exhibited in Figure 5 and Figure 8 for the fitted minus linear theory fits is consistent with 
previous estimates of the reconstruction errors for IRAS type gravity maps (Strauss et al. 
1992). The flat behavior of '^ ^ is indicative of an error in the overall dipole amplitude of 
the IRAS-hke flt, of order 100 km/s, 

It is instructive to examine the value of x^ versus the number of modes included in 
the flt, since this allows an assessment of the statistically signiflcant modes used in the 
reconstruction. The orthogonal function method is very useful for this purpose. The flt 
parameters a^ can easily be arranged in descending order, a^^ > a-'^ • • • > a^rnax ^ g^^^j ^^2 
can be computed as a function of the limiting mode. Figure 9 shows x^ versus the number 
of modes incorporated in the reconstruction, for the case {rimax^ ^max) = (8, 3) in the mock 
Marks catalog constructed with cr^ = 0.03. We see that for a number of modes higher 
than j = 100 the gain in the x^ per degree of freedom is marginal. The "saturation" of 
X^ as a function of degrees of freedom is the criterion according to which the statistically 
signiflcant flt parameters are determined in realistic applications. Note that, due to the 
small scale velocity dispersion which is not fltted by the model, the x^ of the model per 
degree of freedom is larger than unity, even when all the modes are included. 

4. DISCUSSION 

We have presented a method for estimating the peculiar velocities of disklike galaxies 
given measurements of their apparent magnitudes and circular velocity parameters. The 
method depends on a linear approximation of the inverse TuUy-Fisher relationship and 
is free of Malmquist bias in the absence of multi-streaming. The bias induced by the 
observed amplitude of small scale peculiar velocities is essentially negligible. In any case 
the bias can be completely avoided by discarding nearby galaxies in the estimation of 
peculiar velocities. The peculiar velocities are given by means of an expansion in a set of 
independent and orthogonal base functions. Therefore the method successfully recovers 
the smoothed components of the peculiar velocities of galaxies without the necessity of 
prior binning of the data. The choice of the base functions determines the form of the 
smoothing employed. For example, in the description of a discrete set of radial Bessel and 
angular spherical harmonics, which correspond to an expansion of a 3D Fourier series, the 
smoothing is equivalent to a sharp cutoff in A;— space. A general window function may 
be introduced by multiplying the Bessel and spherical harmonics expansion by its inverse. 
The model peculiar velocity fleld can be used to generate smoothed maps of the peculiar 
velocity fleld in regions of space where data is not available, although the inversion of the 
matrix of the orthogonal functions is singular and will require special care. The choice 
of the window function clearly depends on the goal of the application of the method. 
For instance, in order to compare the inverse TuUy-Fisher estimated velocities with the 
IRAS predictions, a variable smoothing is appropriate. On the other hand if we seek 
peculiar velocity maps for comparisons with the predictions of cosmological scenarios, a 
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large constant smoothing is desired. 

The asymptotic flatness of the x^ curve (Figure 9) at a level above that expected 
for a "good" flt to the data is one indication that our particular choice of basis functions 
is not ideally matched to the problem at hand; too many oscillations of the spherical 
Bessel functions occur at large distance where the constraints are few. A larger number of 
basis functions only slightly increases the parameters available for better fltting the local 
flow, while introducing still more unconstrained oscillations at larger distances. A better 
optimized set of basis functions should allow us to reach the same or better x^ value with 
fewer flt parameters. 

The method presented here is designed to give a best estimate for the peculiar veloc- 
ities of galaxies by means of their redshifts. Therefore the associated smoothed peculiar 
velocity fleld is given in redshift space rather than real space, which is ideal for comparison 
to the peculiar velocity fleld computed by Nusser and Davis (1994), which was also in 
redshift space. A somewhat cumbersome property of flows in redshift space is that the 
velocity fleld is not vorticity free when the flow is non-linear. This poses some complica- 
tion when seeking a recovery of the mass-density fluctuations, from the velocity fleld (e.g. 
Nusser et al. 1991). Therefore it is worthwhile to obtain an estimate for the velocity fleld 
in real space. The method proposed here can easily be modifled to iteratively solve for 
the velocity fleld in real space. In a separate paper we intend to test the recovery of the 
velocity fleld in real as well as redshift space using various smoothing window functions 
which are of interest in many applications of the method. The comparison of the smoothed 
velocities derived from the actual MarkS data with the 1.2Jy IRAS predicted velocities is 
underway and shall be reported in due course. 
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FIGURE CAPTIONS 

Figure 1 : The estimated peculiar velocities using the inverse TuUy-Fisher relationship 
vs. the true peculiar velocities for galaxies in the mock IRAS catalog. Only velocities 
of "galaxies" within z = 5000km/s are plotted. The left panel shows the reconstruc- 
tion for (nmaxT^max) = (12,4) while the right panel is for the reconstruction with 
{iT-maxT^max) = (8,3). Velocities are in units of km/s. 

Figure 2 : The difference between the estimated and true peculiar velocities vs. redshift 
of the mock IRAS "galaxies" in figure 1. Redshift is in units of km/s. 

Figure 3 : The estimated vs. true peculiar velocities of "galaxies" in the mock MarkS 
data for three values of cr^ (upper and lower right panels). Also shown (lower left panel) 
the estimated vs. peculiar velocities of the mock MarkS "galaxies" vs. the prediction 
of peculiar velocities from the density fluctuations in the mock IRAS survey. Only 
velocities of "galaxies" within z = 5000 km/s are plotted. 

Figure 4 : The difference between the estimated and N-body peculiar velocities from 
flgure 3 vs. redshift. 

Figure 5 : The ratio of the velocity differences from flgure 4 to the la error in the 
estimated velocities vs. redshift. The errors were computed using equation (13). 

Figure 6 : The estimated and true peculiar velocities in the mock Mark3 "galaxies" on a 
spherical shell between z = 3500 and 5000 km/s. Open circles denote galaxies infalling 
to us (LG frame); stars are outflowing from us. The size of the symbol is proportional 
to the amplitude of the flow; a flow velocity of 200 km/s is indicated in the lower right. 

Figure 7 : The difference between the estimated and true peculiar velocities vs. redshift 
of the mock IRAS "galaxies" . In addition to their true peculiar velocities, the "galaxies" 
were assigned random radial peculiar velocities with dispersion of cr^ = 400 km/s. 

Figure 8 : The square root of the correlation function of velocity residuals (units (km/s) 
versus pair separation. Zero lag corresponds to the rms amplitude of the residuals. 
The solid curve is the correlation of the fltted minus true N-body velocity, while the 
long dashed curve is the correlation function of the fltted minus linear theory velocity. 
The covariance of the true velocity fleld is shown as the short dashed curve. 

Figure 9 : The x^ vs. the number of modes incorporated in the reconstruction of the 
model. The flt parameters were flrst put in descending order. 
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